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0^ I Within the framework of the internal shock model, we study the luminosity 

I and the variability in gamma-ray bursts from collimated fireballs. In particular 

^ I we pay attention to the role of the photosphere due to pairs produced by the 

O ! internal shock synchrotron photons. It is shown that the observed Cepheid-like 

Q I relationship between the luminosity and the variability can be interpreted as a 

O I correlation between the opening angle of the fireball jet and the mass included 

^ I at the explosion with a standard energy output. We also show that such a 

O I correlation can be a natural consequence of the coUapsar model. Narrow jets, in 

'-^1 which the typical Lorentz factors are higher than in wide jets, can produce more 

^ I variable temporal profiles due to smaller angular spreading time scales at the 

I photosphere radius. Using a multiple-shell model, we numerically calculate the 

^ I temporal profiles of gamma-ray bursts and show that our simulations reproduce 

I the observed correlation. 

H ' Subject headings: gamma rays: bursts; shock waves 



1. Introduction 

Gamma-ray bursts (GRBs) and the afterglows are well described by the fireball model 
(e.g. see Piran 2000), in which an explosive flow of relativistic matter (ejecta) is released 
from a central engine. The collision of fast-moving ejecta with slower ejecta results in a 
GRB. Thereafter, the ejecta shocks and sweeps up a large amount of ambient matter. The 
shocked matter powers the long-lived afterglow. 

The temporal profiles of GRBs are often very variable and each profile looks very differ- 
ent. Quantitative measures of the variability have been suggested allowing for a systematic 
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study of their morphology. Several studies have explored the possibility that quantities di- 
rectly measurable in GRB light curves could be related to the luminosities of the bursts. 
Stern, Poutanen & Svensson (1999) concluded that there is an intrinsic correlation between 
luminosity and the complexity of GRBs. Ramirez-Ruiz & Fenimore (1999; see also Fenimore 
& Ramirez-Ruiz 2000) found that the luminosities of seven bursts with known redshifts are 
correlated with the variabilities. Based on this work, Reichart et al. (2001) also reported 
a possible Cepheid-like luminosity estimator for long bursts. Using the GRB rate density 
derived from this correlation, Lloyd-Ronning, Fryer & Ramirez-Ruiz (2001) discussed the 
star formation rate at high redshift. 

Recent afterglow observations allow us to determine the geometry of ejecta, whether it 
is spherical or conical. There is excellent observational evidence that the ejecta have conical 
geometry (jet), and there appears to be a correlation in the sense that the bursts with the 
largest gamma-ray fluences have the narrowest opening angles. The correlation is improved 
when the fluences are all scaled to the same distance by using the redshift measurements. 
Frail et al. (2001), Panaitescu & Kumar (2002) and Piran et al.(2001) reported that the 
gamma-ray energy releases, corrected for geometry, are narrowly clustered around lO^^erg, 
and suggested that the wide variation in fluences and luminosities of GRBs is due entirely 
to a distribution of the opening angles. If this is the case, the Cepheid-like relation implies 
a correlation between the opening angle and the variability: more highly coUimated GRBs 
are more variable. Such a correlation is actually found in the observations (see figure 1). 

As a result of rclativistic beaming, an observer can see only a limited portion of the 
ejecta. There should be no observable distinction between a spherical ejecta and a conical 
ejecta until the ejecta has slowed down in the afterglow phase. Therefore, the correlation 
between variability and opening angle should be attributed to the properties of the central 
engine. In this paper, we will explore the relations among variability, luminosity and jet 
opening angle in the framework of the internal shock model. We show that the variability- 
opening angle relation and the Cepheid-like relation can be interpreted as a correlation 
between the opening angle and the mass included at the explosion. In §2 we first discuss 
the collapsar model as the central engine of GRBs. In §3 we make some comments on the 
typical GRB energy. In §4 we discuss the internal shocks and the temporal profile. In §5 we 
study internal shocks by using a multiple-shell model. In §6 we consider what assumptions 
are actually require to fit the Cepheid-like relation. In §7 we give conclusions. 
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2. Central Engine: Collapsar 

While the nature of the GRB progenitors is still unsettled, it now appears likely that at 
least some bursts originate in explosions of very massive stars (main sequence mass Mms > 
25 Mq). The model currently favoured for long bursts is the collapsar model in which GRBs 
are caused by relativistic jets expelled along the rotation axes of collapsing massive stars. 
The jets are powered by a black hole with a surrounding accreting torus. Energy from the 
accretion is pumped into jets via electrodynamic processes or by neutrino annihilation. In 
addition spin energy of the rotating black hole may be tapped by magnetic fields anchored 
in the accretion disk. 

According to Frail et al. (2000), Panaitescu & Kumar (2002) the wide variation in the 
fluences of GRBs originates from the differences in opening angles of the jets. A wide jet 
radiates gamma-ray photons into a large solid angle, resulting in a dim burst. Though we do 
not know yet what physical mechanism results in the wide variation of the opening angles, 
it is likely that a wider jet involves a larger mass at the explosion, and that it results in a 
flow with a lower Lorentz factor. 

In the collapsar model, the duration of a GRB is set by the collapse time scale of a 
massive stellar core, typically ~ 10 sec for helium cores of 8 < < 15. The progenitor has 
a mass > 25 Mq on the main sequence but loses its envelope to a binary companion and/or 
stellar wind before core collapse. The stellar core is assumed to collapse to a black hole, due 
to it's large mass, and to be rotating sufficiently rapidly to form an accretion disk around 
the newly born black hole (Woosley 1993; MacFadyen & Woosley 1999). 

Typical accretion disk time scales are milliseconds and gas typically resides briefly in 
the accretion disk. The key point is that the disk is continuously fed by the collapsing star. 
Short time scales are available and indeed expected due both to variability in accretion rate 
and the jet instabilities as the hot jet material expands along the polar axis of the star and 
interacts with the surrounding star. 

For disks forming with radius ~ 10 cm, typical accretion time scales are t 

acc ^ 0.01 

sec. Fluctuations in accretion rate due to instabihties were shown to exist with 50 msec 
time scales (MacFadyen & Woosley 1999). In models where accretion energy is tapped to 
power relativistic polar jets, these fluctuations in accretion may translate into variation in 
jet Lorentz factor perhaps leading to variation in the GRB light curves. However time scales 
calculated so far in numerical simulations are probably too short to produce the observed 
variations. 

A more promising source of time structure in GRBs are instabilities arising as a rel- 
ativistic jet propagates through the stellar mantle. Instabilities in the flow due to shear 
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between the jet and the star or backflow from the jet head result in fluctuating jet speeds. 
Recent numerical simulations of relativistic jets propagation through collapsars demonstrate 
the presence of the variability some with time scales of ~ 0.1 sec (Zhang, Woosley & Mac- 
Fadyen) . 

MacFadyen, Woosley & Heger, 2001 (hereafter MWH) showed that jets of equivalent 
energy injected into a stellar envelope can be focussed by the pressure of the star. The 
degree of focusing is a function of the assumed entropy of the jet which was parametrized 
by Eint/Efot- "Colder" jets were squeezed into tightly collimatcd flow while "hotter" jets of 
the same energy expanded sideways (Sec MWH Fig 10). The rcsTilt was that "hotter" jets 
are broader and sweep up a larger mass of stellar material along the rotation axis of the 
star. While the calculations in MWH were non-relativistic, the results should hold in special 
relativity since the asymptotic Lorentz factor is an inverse function of the swept up mass. 
Current fully relativistic numerical simulations indicate the results do hold (Zhang, Woosley 
& MacFadyen). 

MWH also experimented with the effect of the coUapsar density structure on jet coUima- 
tion and found a large difference in jet opening angle for two density structures corresponding 
to low and high disk mass. The difference was due to the assumed viscosity parameter, but 
may also result from initial angular momentum and density structure of the progenitor stars. 
Calculations of the collapsar density structures and their effect on relativistic jet opening 
angle are currently underway. 

3. Typical GRB energy 

In the collapsar model, the energy powering the GRB originates in an accretion torus 
surrounding a rotating black hole of several solar masses. The stellar progenitor must lose its 
envelope in order to have a sufficiently small radius to allow a relativistic jet formed near the 
accreting black hole to pierce the star in a GRB time scale i.e. R^, < lOsec (c/2) ~ 1.5 x 10^^ 
cm. In addition the range of angular momentum conducive to GRB formation may be narrow. 
It is constrained on the low end by the angular momentum sufficient to form a disk and on the 
upper end, so that the disk is small enough to cool at least partially by neutrino emission 
(though these neutrino need not participate directly in jet formation). Disks which are 
inefficiently cooled are convective leading to outflows and inefficient accretion (MacFadyen 
& Woosley 1999; Narayan et. al). The successful GRB-producing collapsars may therefore 
be expected to have roughly similar energetics. 

Another requirement for GRB production is that the jets have sufficient momentum flux 
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to overcome the ram pressure of the column of gas accreting along the stellar rotation axis. 
Since the momentum flux decreases with a larger opening angle, there is a maximum opening 
angle which will not be swamped by the accretion ram. Equivalently for a given opening 
angle there is a minimum momentum flux capable of launching a jet. Given an observed 
maximum jet opening angle, this sets a lower limit for energy flux capable of producing a 
GRB. Given a similar typical time scale this implies a similar total energy. 

4. Internal Shocks and the Temporal Structure 

Internal shocks arise in a relativistic wind with a nonuniform velocity when the fast 
moving flow catches up the slower one. The wind can be modelled by a succession of 
relativistic shells (Kobayashi, Piran & Sari 1997, hereafter KPS). A collision of two shells 
is the elementary process, and produces a single pulse of gamma-ray. Three time scales, 
the cooling time, the hydrodynamic time and the angular spreading time, are relevant to 
the pulse width. With the relevant parameters the cooling time is negligible compared to 
the other two time scales (Sari, Narayan & Piran 1996). Let d and D be the width and 
separation of the shells. In other words, we assume that the central engine operates for a 
period d/c and then is quiet for a period D/c. The hydrodynamic time scale d/c and 
the angular spreading time scale ^ Djc determine the rise and the decay time of the pulse, 
respectively. Since most observed pulses rise more quickly than they decay, the pulse width 
is mainly determined by the angular spreading time ^ D/c (Norris et al. 1996 ; Ryde & 
Petrosian 2001). 

The whole light curve of a GRB is given by the superposition of the resulting pulses 
from each collision. Since all shells are moving towards us with almost speed of light, we 
observe pulses arising from collision between shells mostly according to their positions inside 
the wind, i.e., according to the time when those shells were emitted by the central engine 
(KPS; Nakar & Piran 2001). The relative positions of the shells inside the wind are also 
determined by the separations. Therefore, the variability time scales in the temporal profile 
reflect the shell separations at the central engine well. 

The observed distribution of peak separations {D/c} has a large dispersion (Nakar & 
Piran 2001). While the fastest variability time scale in GRBs is about a miUisecond, the 
largest peak separations or the total durations are usually very much longer. Long bursts 
have durations from 10 to lOOOsec. Clearly, whatever is the physical mechanism behind the 
GRB production, it acts on a much longer timcscalc than the fastest dynamical time. Many 
GRB models involve accretion onto a compact object, usually a black hole. If the accretion 
disk is fed by fallback of material after a supernova explosion, as in the coUapsar model. 
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then the duration and large peak separations are determined by fallback, not the dynamical 
accretion time. The large dispersion of peak separations could be explained in the coUapsar 
model. 

It has recently been shown that the Thomson optical depth due to pairs produced 
by synchrotron photons plays an important role in the internal shock model (Guetta, Spada 
& Waxman 2001; Asano & Kobayashi 2002). In order to obtain high radiative efficiency and 
the characteristic clustering of spectral break energies of GRBs in the range 0.1-1 MeV, the 
collision radii arc required to be similar to the photosphere radius. Since collisions producing 
narrow pulses occur at small radii oc D, the photosphere might obscure these, leaving only 
the wider pulses visible. This will make the temporal profile smooth. In a wider jet, the 
typical Lorentz factor F is smaller, a larger fraction of the collisions occur at small radii 
oc below the photosphere, therefore, the smoothing effect is expected to be stronger. This 
might explain the correlation between luminosity and variability in GRBs. 



We discuss the smoothing effect by using a multiple-shell model. We represent the 
irregular wind by relativistic shells in a manner similar to that in KPS. Because of the 
relativistic beaming effect, we can study the emission from a jet by using spherical shells 
with an isotropic explosion energy Eiso{0) — ^E/O"^ where E is the geometrically corrected 
explosion energy and 9 is the opening angle. The gamma-ray energy released from a GRB 
is narrowly clustered around lO^^erg (Frail et al. 2001), and the conversion efficiency from 
the explosion energy into the gamma-rays is about 10% (Guetta, Spada & Waxman 2001; 
Kobayashi & Sari 2001). Then, E ~ 10^^ gj,g 



Consider a two shell collision which is the elementary process in the multiple shell 
evolution. A rapid shell with Lorentz factor and mass catches up to a slower one 
with F^ and m^, and the two merge to temporarily form a single shell. Using conservation 
of energy and momentum we calculate the Lorentz factor of the merged shell to be 



The internal energy of the merged shell is the difference of kinetic energy before and after the 
collision, E^nt — 'mrC^{Vr—^m)-\-'^sC^iys—^m)- We assumc that electrons are accelerated to a 



5. 



Shell Model 



5.1. Two-Shell Collision and Photosphere 




(1) 
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power-law distribution of Lorentz factor 7e, with a minimum Lorentz factor 7^: N{'^e)dle 
7~^(i7e,7e > 7m- Throughout this paper, we use the standard choice p = 2.2. A constant 
fraction of the internal energy goes into magnetic energy. The energy distributed to the 
electrons is radiated via synchrotron emission. 

A significant fraction of the kinetic energy could be converted to the internal energy if 
the relative velocity of the two shells is relativistic. However, it is not reasonable that all 
the internal energy is emitted, because electrons do not have most of the internal energy, 
and protons do not cool. Defining eg as the fraction of energy given to electrons, we expect 
ee < 1. Even at equipartition with protons eg is 1/2 

If ee < 1, the merger produced by a collision is expected to stay hot after the emission. 
As a result, the merger will spread, transforming the remaining internal energy back to 
kinetic energy. Kobayashi and Sari (2001) studied such a process by using a hydrodynamic 
code to show that the spreading leads to a formation of two shell structures, rather than a 
homogeneous wide shell. This is a planar analogy of the evolution of a relativistic fireball 
(Kobayashi, Piran and Sari 1999). 

Let us consider a two shell collision in the rapid shell rest frame. When the slow shell 
begins to interact with the rapid shell material, two shocks are formed: a forward shock 
propagating into the rapid shell and a reverse shock propagating into the slow shell. After the 
reverse shock crosses the slow shell, the profile of the shocked rapid shell material approaches 
to that of its fireball analogy: the "blast wave" which sweeps and collects the rapid shell 
material. Once the shock wave crosses the rapid shell, the hydrodynamical structure is as 
follows. There is a shocked rapid shell, the analogy of the "blast wave" and the shocked slow 
shell, which is the analogy of the "fireball ejecta" . 

When the evolution of multiple shells is consider in the following sections, shell collisions 

happen around a radius ~ T^D, and therefore the time interval of the collisions is order of 
~ r'^D/c. Since the hydrodynamic time scale (shock crossing time) ~ T^rf/c is much shorter 
than the interval, each two shell interaction is completed in a moment, and no shell collides 
into the merger during the interaction. A simplified description of the spreading effect is to 
assume that the two shells reflects with a smaller relative velocity (Kobayashi & Sari 2001). 

The synchrotron photons can be scattered by electrons within the merger. The Thomson 
optical depth is increased significantly when taking into account pairs produced by internal 
shocks (Guctta, Spada & Waxman 2001; Asano & Kobayashi 2001). The typical energy of 
the synchrotron emission, in the shell frame, is 

hu'^ ^ 10 kev ey^'Elil,,Ri,%'/'r;;,],, (2) 

where / = 10^/^;, is in units of erg, the collision radius R and the shell width d are in 
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units of cm. This energy is well below the pair production threshold. However, since the 
photon spectrum extends to high energy as a power law 

there exists a large number of photons beyond the threshold hv' ~ mgC^. The pairs produced 
by these photons can contribute significantly to the Thomson optical depth. In order to take 
the effect of pair production into account, we determine for each collision the typical photon 
energy hv'^ and the amount of energy emitted above the threshold. 

, , X -p/2+1 p 



E' ~ 1^ ^"^"^ (4) 

Since the number density of the pairs is given by n'^ ~ E'^/AnnieC^ R^dTm, the photosphere 
radius R±, where the Thomson optical depth becomes unity, is given by 

If a collision happens below the photosphere, the whole internal energy produced by the 
collision is converted to kinetic energy again via the shell spreading. We assume that the 
two shells just reflect with the same relative velocity. 



5.2. Four Shell Evolution: an example case 

We here consider a evolution of four shells which is a basic component in multiple shell 
interaction. Since the observations require that shell widths should be much smaller than 
the separations, we assume infinitesimally thin shells for simplicity. We assign an index i, 
(i = 1,4) to each shell according to the order of the emission from the central engine (see 
figure 2). The central engine ejects each shell with a Lorentz factor Fj 3> 1 at time Tj. We 
choose the origin of time such that the first shell is emitted at Ti = 0. For convenience we 
define A = c(Tj — Ti_i) as the separation between shell i — 1 and shell i. 

To evaluate the conversion efficiency or the temporal properties of the internal shock 
process, it is important to understand how high velocity shells interact with slower ones and 
dissipate kinetic energy. Thus, we consider a case that the Lorentz factors of the shell 2 and 
the shell 4 are much higher than those of the shell 1 and the shell 3 F2, F4 » Fi, F3. The 
collision between the shell 1 and the shell 2 takes place at time T ~ 72 + 2FfD2/c and radius 
R ~ 2FfD2. The gamma-ray pulse produced by this coUision arrives at the observer at the 
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observer time t^T- R/c ~ + o(rf/r|). We chose the origin of the observer time such 
that if a pulse was emitted from the engine at T = 0, it would arrive at the observer at t = 
(dashed dotted line in figure 2). The pulse produced by the collision between the shell 3 and 
the shell 4 comes to the observer at t ~ T4 + o(r|/r|). 

The mergers (filled circles) spread and cause a subsequent collision (open circle) . Using 
the Lorentz factor of the shell 2 and the shell 3 after the collisions and Fg, the arrival 
time of the pulse from the subsequent collision can be estimated as t ~ T4 + o{T'^/T'^). If 
F2/F3 ^ 1, the pulses from the collisions between the shell 3 and the shell 4 and between 
the shell 2 and the shell 3 arc observed almost at the same time. If F2 ~ F3, the collision 
can produce only a weak and soft pulse, it does not contribute to the main feature of the 
GRB light curve. Therefore, we observe pulses arising from the collisions mostly according 
to the time when high velocity shells were emitted by the central engine. 



5.3. Multiple-Shells 

We consider a wind consisting of shells. Each shell is characterised by four variables: 
Lorentz factor Fj, mass rrii, width di and the distance to the outer neighbour shell Di. As we 
wrote in section 4, most observed pulses in GRB light curves rise more quickly than decay, 
this means that the widths are smaller than the separations. Since, in such a case, light 
curves are not sensitive to the distribution of widths, for simplicity the width is assumed to 
be constant d = Dmin- Though we will show numerical results with various distribiitions 
of Fj, rrii and later, we first consider simple distributions with which we can give simple 
arguments to understand the numerical results. We assume that the Lorentz factors and 
the separations are distributed uniformly in logarithmic spaces; between F„jj„ and r^ax 
and between Dmin and Dmax, respectively. We assume that the shells have equal mass. The 
internal shock process becomes very efficient if shells have the same mass and if the dispersion 
of the initial Lorentz factors is large (Beloborodov 2000; Kobayashi and Sari 2001). The value 
of mass is normalised by using the isotropic explosion energy Eigo{9) sts m — E'j^o/SFjC^. 

Our key assumption is a correlation between the Lorentz factors and the opening angle: 
wider jets have lower Lorentz factors F oc 9"°'. The observed wide jets with 9 ~ 0.2 also 
should have ultra-relativistic Lorentz factors to be optically thin to high energy photons, we 
assume r^axiO) = ^0 (9/0.2)-'^ and Fq = 100 

The evolution of the system in time is basically the same as in KPS97, but the photo- 
spheric effect is considered. We follow the evolution of shells until there are no more collisions, 
i.e. until the shells are ordered by increasing values of the Lorentz factors. The numerical 
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light curves in the 50 - 300 keV BATSE band are plotted in Figure 3a-c for a model with dif- 
ferent opening angles 6 = 0.2, 0.06 or 0.02. Since we have assumed a = 1, these correspond 
to T^ax = 100,333 and 1000, respectively. = 100, Tmax/^min = 10, -Dmm/c = 1msec, 
Dmax/c = Isec and eg = = 0.1 were also assumed. As we expected, the temporal profiles 
are more variable for narrower jets with higher Lorentz factors. 

The definition of variability by Fcnimore & Ramirez-Ruiz (2000) and that by Reichart 
et al (2001) are slightly different, but both relate the variability to the square of the time 
history after removing low frequencies by smoothing. To evaluate the variability of numerical 
light curves, we use a simplified version of the Reichart et al. (2001) variability, 

E(Q-(C))! 

where Cj is 50-300 keV photon counts at time ti with 64 msec resolution (the BATSE 
resolution), and (C) is the counts smoothed with a boxcar window with a time scale equal 
to the smallest fraction of the burst time history that contains a fraction / of the total 
counts. Reichart et al. (2001) found that / = 45% gives a robust definition of variability. 
Using the formula (6) with / = 0.45, we evaluated the variability measures asV — 0.06, 0.10 
and 0.25 for Figure 3a, 3b and 3c, respectively. 

Reichart et al (2001) found that the isotropic peak luminosities L are correlated with 
the variability measures, L oc We now show that such a correlation exists in our 

numerical simulations also. A numerical temporal profile depends on the specific realizations 
of random variables that are assigned to each shell, e.g. random Lorentz factors and random 
separations, as well as on the model parameters. For a given 6 and the same model param- 
eters as in Figure 3, we calculate the temporal profiles for 100 realizations, and evaluate the 
mean isotropic peak luminosity and the mean variability measure. Here the peak luminosity 
is calculated with 1 sec resolution as Reichart et al (2001) had analysed. In Figure 4a, the 
thick solid line shows the mean values, while the thin solid hues depict the la error. We also 
plot the opening angle- variability relation in Figure 4b. The numerical results reasonably fit 
the observational data, and can be understood using the following arguments. 



5.4. Analytic Estimate 

Numerous collisions happen during an evolution of multiple-shells. Each collision pro- 
duces a pulse. However, the main pulses are produced by collisions between the fastest 
shells and the slowest shells ~ l^min- Since such a coUision happens at it! ~ V^^^Di 

for an initial separation Dj, we can define two characteristic collision radii, F^-^D^in and 
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T'^^^Dmax- In this equal mass case, mass is given by mc^ ~ 2Eiso/NTmax- The internal 
energy produced by a main collision and the Lorentz factor of the merger are Ei^t ~ Eiso/N 
and ~ \/^max^mim respectively. The Thomson photosphere radius due to pairs from 
pair production interaction of synchrotron photons can be estimated by assuming that a 
significant fraction ~ 10% of the radiative energy ^eEint is converted to pairs. 



~ 5 X 10^^ cm Eli'N-'/'T,^J\e/0.2)-'/\ (8) 
hereafter eg = 0.1 is assumed. 

If R± > r'^-^Djnax: the first coUisions for most shells occur below the photosphere. How- 
ever, since the shells reflect each other with the same velocity, it is still possible to produce 
bright pulses when the shells propagate beyond the photosphere and collide into other shells. 
If an enormous number of collisions and reflections happen below the photosphere, the shells 
are ordered with increasing values of the Lorentz factors and no collision happens any more. 
In order to produce bright pulses, the shells should come out from the photosphere before 
the slowest shells are sorted out at the tail of the wind with a total width A ~ ND„^ax/2. 

The equalities, R± = V^^^Dmin and R± = F^^^A give two opening angles. 

0, ~ 10-' E;,'/'N',/'rl/^,'D'X_, (9) 

e, ~ 10-' E;,'/'N,rl/^D'±^, (lo) 

where D„un,~3 = (/^mm/c) /lO^^sec and Dmax,o = {Dmax/c) /Isec. Using these opening 
angles, GRBs are classified into the following three cases. (1) Narrow Jet Case: 6 < 9i. 
Since all collisions occur above the photosphere, the variability should be less dependent on 
6 and only on the distribution of the separations at the central engine. (2) Intermediate Jet 
Case: 9i < 9 < 62. The variabihty of the temporal profile should highly depend on 9. Some 
of the main collisions occur above the photosphere to produce bright pulses wider than 

5t ^ i?±/2crL„ - 10-^sec Eli'N-'^'T-;f{9/10-'f/\ (11) 

(3) Wide Jet Case: 9 > 92- The slowest shells are sorted out at the tail of the wind within 
the photosphere, only minor collisions happen. The resulting bursts become dimer. 

Figure 3a is an example of the Wide Jet Case, while figure 3b and 3c correspond to the 
Intermediate Jet Case. Equation (11) gives rough estimates on the widths of main pulses for 
these Intermediate Jet Cases, 5i ~ 3 sec {9 — 0.06) and ~ 0.5 sec {9 — 0.02). Main pulses are 
actually wider in figure 3b. Prom figure 4b, the variability (thick solid line) is about 0.08 at 
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— 02 0.1. In figure 4a, a linear fit to the numerical data for V > 0.08 is L oc V^-^, which 
is consistent with a power law reported by Fenimore&Ramirez-Ruiz(2000) and by Reichart 
et al. (2001). Since the shells are ordered by velocity within the photosphere, the numerical 
luminosity rapidly decreases for V < 0.08. 

6. Other Initial Conditions 

In our numerical simulation, we have assumed that a wider jet has a lower Lorentz 
factor r oc 0~°'. A second assumption is that the shells have equal mass. A third concerns 
the distributions of the initial separations and the initial Lorentz factors. In this section, we 
discuss how numerical results depend on these assumptions. 

Since the dynamics of the blast wave is determined only by two parameters: isotropic 
explosion energy and ambient matter density, afterglow emission does not depend on the 
initial Lorentz factor. The current afterglow observations detect radiation from several hours 
after the bursts. At this stage, the Lorentz factor of the blast wave is less than ~ 10. 
Then, the initial Lorentz factor can not be directly constrained through afterglow modelling. 
However, the Lorentz factor at the beginning of the afterglow stage can be estimated if one 
knows when the afterglow began. 

An optical flash was observed during an exceptionally bright burst GRB990123, by 
using this peak time, the Lorentz factor was estimated as about 300 (Sari & Piran 1999; 
Kobayashi 2000). However, the optical flash was detected only for it so far. Panaitescu 
and Kumar (2002) assumed that the observed GRB duration is a good measure of the peak 
time, and found a correlation F oc within ten events for which opening angles are well 
constrained from the afterglow break time. Salmonson and Galama (2002) also have recently 
suggested the anticorrelation of F with 0, based on the positive correlation they observed in 
several cases between the GRB pulse lag time and the afterglow break time. However, the 
dependence they infer is much stronger F oc 0~^/^. Therefore, we have done simulations with 
a = 1/2 and 2 also. 

The numerical results with a = 2 are shown in figure 4a and 4b (dashed lines). The 
other parameters are the same with the case of a = 1 (solid line). For a larger value of 
«, the Lorentz factor depends more strongly on the opening angle, while the photosphere 
radius becomes insensitive R± oc ^("~2)/2 gj^ce the variability V increases more rapidly 
with decreasing of the opening angle, the variability - luminosity relation becomes gentle. 
For a = 1/2, the numerical variability is insensitive to the opening angle. When we calculate 
y in a range of 0.02 < < 0.5, it is almost constant V ~ 0.05. Therefore, the variabihty - 
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luminosity relation becomes very steep. The numerical results with a ~ 1 give the best fit 
to the observations. The mass loading at the central engine is M ~ E/c^T oc 6, which is in 
a good agreement with a correlation M oc 9^'^ reported by Panaitescu and Kumar (2002). 

We next discuss the assumption about mass. We also studied the cases that the shells 
have the same energy initially or have random masses. The random masses are assumed to 
be uniformly distributed between m and 10m. The value of m is normalised at the beginning 
of the computation by the isotropic energy Eiso = TiiTLiC^Vi. Under these mass assumptions, 
we calculate the temporal profiles for 100 realizations, and evaluate the mean isotropic peak 
luminosity and the mean variability measure. Figure 5a and 5b present the variability - 
luminosity and variability - opening angle relations. The equal energy (dashed line) and 
random mass (dashed dotted line) cases give similar results to the equal mass case (solid 
line). For a given 9, the variability in the random mass case is larger than in the equal 
mass case. However, by choosing a smaller normalization of Lorentz factor Fq ~ 80, we can 
shift the line horizontally, and can get a more similar relation. Therefore, the equal mass 
assumption is not essential to get the observed correlation. 

GRB light curves contain many pulses with widths of msec to sec time scales. Since 
widths are determined by the initial separations, we considered a simple case that the initial 

separations Di are uniformly distributed in logarithmic spaces between 1 msec and Isec. 
However, it is known that pulse widths have a lognormal distribution (Nakar & Piran 2001). 
We here consider a lognormal random intervals: a mean (log^Q(Z})) = — 2 and a dispersion 
5{\og^Q{D)) = 1. The thick solid line in figure 6 shows the numerical results, and is similar 
to the case of uniform distribution (thin solid line). 

In order to achieve a high conversion efficiency (> a few%) from the explosion energy to 
gamma-ray in the internal shock process Tmaxf^min should be larger than ~ 10 (Kobayashi 
& Sari 2001). The dashed dotted line in figure 6 depicts the case of Tmaxf^min = 100, which 
steepness is comparable to the case of Tmaxf^min = 10 (thin solid). For a given luminosity, 
the variability is larger in this case. However, by using a smaller Fq, we can get a better fit to 
the observations. We also have done numerical simulations with a Gaussian random Lorentz 
factors where a mean value (F) = 50(^/0.2)"^ and a dispersion ST — (F) are assumed. We 
truncated the initial Lorentz factor distribution for F < 2, e.g. when we assign a random 
value to each shell at the beginning, random Lorentz factors are generated until F, > 2 is 
satisfied. The numerical results (dashed line) are also similar to the other cases (see figure 
6). 

We have shown that the variability - luminosity relation exists even if we make different 
assumptions on the distributions of masses, intervals and Lorentz factors and that the index 
and the normalization of the variability - luminosity relation highly depends on a and Fq. 
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Our numerical results suggest a ~ 1 and Fq ~ 10^. 

Another assumption concerns the values of and e^, fractions of internal energy dis- 
tributed to electrons and magnetic field at coUisions. We have done another set of simulations 
with Ee — €3 — 0.3 and get very similar results. 

7. Conclusions 

We have shown that there exists a correlation between the jet opening angle 9 and the 
gamma-ray light curve variability V. Though the correlation is based on only seven events 
at present and needs to be further confirmed with more events, it is naturally expected if 
the luminosity L is correlated with the variabihty (Fenimore & Ramirez-Ruiz 2000; Reichart 
et al. 2001), and if GRBs have a standard energy output (Frail et al. 2001; Panaitescu 
& Kumar 2002; Piran et al. 2001). This correlation might give us a way to measure the 
opening angle for a long burst directly from the GRB light curve. 

We have shown that the opening angle - variability relation, or equivalently, due to 
the constancy of burst energy, the luminosity - variability relation can be interpreted as the 
correlation between the opening angle of a fireball jet and the Lorentz factor. Larger opening 
angles are expected to be associated with greater mass loading at the central engine and 
might result in lower Lorentz factors. We also show that such a correlation can be a natural 
consequence of the coUapsar model. Using a multiple-shell model, we numerically calculate 
the temporal profiles and estimate the luminosity and the variability. Our numerical results 
suggest F oc ^"■'^ or equivalently M oc 6. 

If the opening angle of a jet is very wide 9 ^ 92, the shells are almost ordered by 
increasing values of the Lorentz factors in the photosphere, only minor collisions or, at the 
extreme, only external shocks happen. The resulting dim burst should be smooth and has a 
soft spectrum. GRB980425 and recently reported Fast X-ray Transients (Heise et al. 2001; 
Kippen et al 2001) could be classified into this "Wide Jet Case". 

Norris, Marani & Bonnell (2000) found that the luminosity L is inversely proportional 
to the GRB pulse lag time. The detailed study on the lag time is beyond the scope of this 
paper. However, if the time lag is also determined by the angular spreading time scale, which 
is the key time scale in our model, we can get the same relation. In our model, luminosity is 
scaled by the opening angle L oc 9'^, while the angular spreading time is tang F~^ oc 9^. 
Then, we obtain L oc 

Another measure of variability in light curves is the power density spectra (PDSs). 
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The observed PDSs show a strong correlation with the luminosity (Beloborodov, Stern and 
Svensson 2000). In each case of figure 3 a, b and c, we calculate the light curves for 100 
realizations, and evaluate the average PDS Pf. The averaging means that we sum up the 
PDSs of the numerical profiles with the peak normalisation, and divided by the number of 
realizations, 100. The numerical PDSs show a strong correlation with the opening angles 
: brighter bursts with smaller opening angles have more variable power at high frequency 
(see figure 7). It resembles the results for different brightness classes of real bursts (figure 
8 in Beloborodov, Stern and Svensson 2000). The numerical results reproduce the slope 
index of ~ —5/3 and the 1 Hz break also which Beloborodov, Stern and Svensson (2000) 
found in observed light curves. The photosphere radius R± is about 5 x 10^^ cm and weakly 
depends on parameters. The collisions which produce Iscc pulses happen around a radius 
~ 6 X 10^^(rmm/100)^ cm. If Tmin ~ 100, the IHz break can be explained by the phtospheric 
effect. The constraints on the activity of the central engine from the PDS study will be 
discussed in the subsequent paper. 

In our model, more variable bursts have higher Lorentz factors, and therefore are ex- 
pected to have higher spectral peak energy. We can actually see such a trend in the ob- 
servational data when correcting the cosmological expansion effect (Table 1 of Fenimore and 
Ramirez- Ruiz 2000). Ramirez-Ruiz and Lloyd- Ronning (2002) recently extended our study 
to examine the spectral property of GRBs. Their numerical results reproduce the trend. 

Most studies on GRBs assume a jet with a well-defined angle, though this angle differs 
for different bursts. However, Rossi, Lazzati & Recs (2001) and Zhang & Mcszaros (2001) 
proposed an alternative model where the jet, rather than having a uniform profile out to 
some definite cone angle, has a beam pattern. Even in this model, the more luminous part of 
a jet has a higher Lorentz factor. Then, our results can be applied to explain the variability 
- luminosity relation. 
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Fig. 1. — Opening Angle vs. Variability: Measurements (filled circles) and upper or lower 
limits (arrows) are shown from Reichart et al (2001) and Frail et al (2001). A linear fit to 
filled circles gives 9 oc V'^-^. 
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Fig. 2. — Emission from collisions. Trajectories: shells ejected from the central engine 
with hight Lorentz factor (thick solid), ones ejected with low Lorentz factor (thin solid) 
and observer (dashed). Collisions (circles) produce pulses of photons (arrows). The dashed 
dotted line shows the null geodesic passing through the central engine at T = 0. 
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Fig. 4. — (a) Peak Luminosity vs. Variability with 1 a error bars of 100 random simulations. 
Observational data: measurements (filled circles) and upper or lower limits (arrows) are 
shown from Reichart et al (2001) and Frail et al (2001). a — l(sohd line), a — 2(dashed 
line) (b) Opening angles vs. Variabihty 
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Fig. 5. — (a) Peak Luminosity vs. Variability : mean of 100 random simulations. Equal 
mass (solid line), equal energy (dashed line) and random mass (dashed dotted hne). The 
parameters are the same with Figure 4(q; = 1) except the mass conditions, (b) Opening 
angles vs. Variability 
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Fig. 6. — (a) Peak Luminosity vs. Variability : mean of 100 random simulations. Lognormal 
random intervals (thick solid). T^ax/^min = 100 (dashed dotted) and Gaussian random 
Lorentz factors (dashed), rjnax/^min — 10, uniform distributions of intervals and Lorentz 
factors are assume if the value or the distributions are not specified, a — 1 and equal mass 
are assumed in all cases. The thin sohd line depicts the standard case. 
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Fig. 7. — The averaged PDS for 100 random realizations, opening angle 9 = 0.2(dashed 
dotted), 0.06 (solid) and 0.02 (dashed). The initial parameters are the same as in Fig 3 



